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CD ' Geometrical properties of holonomic and non holonomic varieties defined by the Pfaff equations 

connected with the first order system of equations are studied. The Riemann extensions of affme 
' connected spaces for investigations of geodesies and asymptotic lines are used 

S: 

1 Introduction 

There is the connection between of the Pfaff equation 

P(x,y,z)-^+Q(x,y,z)^+R(x,y,z)^ = (1) 
and the first order system differential equations in form 

O ' dX = dy = dZ (2) 

' P{x,y,z) Q{x,y,z) R(x,y,z)' 

> 

£h ! dP(x,y,z) dQ(x,y,z) = dQ(x,y,z) dR(x,y,z) = 

dy dx dz dy 



> 

^ ' dR(x,y,z) _ dP(x,y,z) 



dx dz 

and its general integral determines the family of the surfaces in R 3 space 

V(x, y, z) = constant, 
which arc orthogonal to the lines of the vector field 

TV = (P(x, y, z), Q(x, y, z),R(x, y, z)). (3) 

In more general case 
the equation (1) 

fi (P{x, y, z)dx + Q{x, y, z)dy + R{x, y, z)dz) = dU (x, y, z) 



1 INTRODUCTION 



2 



is also integrable by means of integrating multiplier a and determines the family of the surfaces 
U(x,y,z) — const passing through the each point of the space and orthogonal to the vector field 
(3). ' ' 

The vector field (3) in the space R 3 with conditions 



(fi, rotN^j = 



is called holonomic. 

In general case the Pfaff equation (1) is not integrable and them corresponds the system of the 
integral curves (Pfaff variety) passing through each point (x, y, z) with tangent lines lying on the plane 



P(x, y, z)(X -x) + Q(x, y, z)(Y - y) + R(x, y, z)(Z -z) = 0. 



(4) 



The set of the planes and points (4) defined by the equation (1) in general case forms two-dimensional 
non holonomic variety M 2 which is generalization of the surface. 

For the variety M 2 may be extended many of the results of classical differential geometry of surfaces. 
For example the notion of the asymptotic lines, curvature lines and geodesic has analog for the variety 
M 2 [1, 2, 3]. 

In fact the solutions of the system 

Pdx + Qdy + Rdz = 0, dPdx + dQdy + dRdz = 



or 

Px + y + Rz = 0, 

P x {x) 2 + Q v {yf + R z {zf + (P y + Q x )xy + {P z + R x )xz + (Q z + R y )yz = 

give us the curve lines of the variety M 2 which are the analog of the of asymptotic lines on the holonomic 
surface. 

The notion of the curvature lines also can be generalized on the variety M 2 . 

They may be of the two kinds and one of them is defined by the solutions of the system of equations 

Pdx + Qdy + Rdz = 0, 

2P x dx + {Q x + P y )dy + (P x + R x )dz P dx 
(Q x + Py)dx + 2Q v dy+(Q z + Ry)dz Q dy =0. 
(P z + R x )dx+(Q z +Ry)dy + 2R z dz R dz 

The notion of the geodesies also can be extended on the variety V 2 (1) and they may be of two 
types. 

The first type is determined from the condition 



P(x,y,z) Q(x,y,z) R(x,y,z) 



d_ 

ds 



x(s) 



ds 



d: 



As) 



d 2 d 2 d 2 

L ds~ 2X{s) ds^ yis) ds^ Z{s) 



0, 



or 



p (x,y,z)—y(s) 



ds 



x(s) Q(x,y,z] 



ds 2 



z(a)+ 



ds 



+ ( -P(x, y, z)—z(s) + ( —x(s) ) R(x, y, z) ) y(s)+ 



ds 
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This relation is equivalent the system of equations 



d 2 x P fdxdP dydQ dzdR\ 



ds 2 P 2 + Q 2 + R 2 \ ds ds ds ds ds ds 

d 2 y Q / dx dP dydQ dz dR\ 

~cb 2 + P 2 + Q 2 + R 2 \ds~~ds~ + ds~~ds~ + ds~~ds~ ) = ' ( ' 

d 2 z R fdxdP dydQ dzdR\_ 

ds^ + P 2 + Q 2 + R 2 \ds~~ds~ + ds~~ds~ + ds~~ds~ J ~ ' ^ 

Remark that after substitution of the corresponding expressions for the second derivatives on coor- 
dinates from the (5)-(7) into the relations 

dP(x,y,z) dx | dQ{x,y,z) dy | dR(x,y,z) dz | z)—+ 
ds ds ds ds ds ds ' ds 2 

d 2 y d 2 z 

+Q(x,y,z)-^ + R(x,y,z)-j-£ = 

one get the identity. 

The definition of the second type of geodesic in the system of integral curves of the Pfaff equation 
(1) is more complicated and does not be used hereinafter. 

2 The Riemann extension of the affine connected space 

The formulas (5)— (7) can be rewritten in form of geodesic of the i? 3 -space equipped with symmetrical 
affine connection W- k {x l ) = W k -{x l ) 

d 2 x l dx j dx k 

+ n fe -7r-7T = °- 



ds 2 Jk ds ds 
In our case we get the system of equations 



d 2 x P ( ( dx\ 2 , „ s dxdy , ^dxdz ^ f dy^ ' 

+ a [ p - U) +(p " + + s +<3 » U 
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where 

A = P 2 + Q 2 + R 2 , 

from which we get the expressions for the coefficients of affine connection. They are 

PPx PQy -j-fl PRz 



TT 1 — — TT 1 — v TT 

— A i 1X 22 — ^ ' 11 33 ~ ^ i 

n l _ P(Py+Q X ) nl _ P^+A,) j _ 

"12- 2A . "13- 2A . "23- 2A 

2 _ Q-Pg tj2 _ QQy tt2 _ QRz 

n n — ~ 22 — a ' ~~ A ' 

n 2 _ Q(P, + Ox) n2 _ Q(P,+Px) n2 _Q{Qz+Ry) 

12 ~ 2A ' 13 ~ 2A ' 23 ~ 2A : 

3 _ RP X 3 _ RQ y 3 _ RR Z 

11 ll — A ' 22 — A ' 33 



A ' 22 ~ A ' 33 ~ A ' 

n3 _ R(Py + Qx) n3 _ R(P Z + R X ) n3 _ R{Q Z + Ry) 
12 ~ 2A ' 13 ~ 2A ' 23 ~ 2A ' 



So with any equation (1) can be associated 3-dimcnsional affine connected space and its properties 
will be dependent from the coefficients of connection 11^ ■ which are determined by the functions P, Q, R. 

In general case such type of connection is not metrizable and corresponding space is not a Ricmannian. 

Further we apply the notion of the Riemann extension of nonriemannian space which were used 
earlier in [4, 5, 6]. 

Remind basic properties of this construction. 

With help of the coefficients of affine connection of a given n-dimensional space can be introduced 
a 2n-dimensional Riemann space D 2n with local coordinates (x % , having the metric in form 

6 ds 2 = -2U^(x l )^f k dx i dx j + 2dt> k dx k (8) 

where ^ k are the additional coordinates. 

The important property of such type metric is that the geodesic equations of metric (8) decomposes 
into two parts 

x k + T k j x i x j = 0, (9) 

and 

-^± + R l kJi Alx i *, = 0, (10) 

where 

SVk d^k_ i dxi 
ds ds jk 1 ds 

and R k ^ are the curvature tensor of 3-dimensional space with a given affine connection. 

The first part (9) is the system of equations for geodesic of basic space with local coordinates x % and 
they does not contains the supplementary coordinates 

The second part (10) of the full system of geodesies has the form of linear 3x3 matrix system of 
second order ODE's for supplementary coordinates \ffc 

d 2 *l> d^ 

_ + ^)_ + B(,)* = 0. (11) 
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It is important to note that the geometry of extended space connects with geometry of basic space. 
For example the property of the space to be Ricci-fiat keeps also for the extended space. 

This fact give us the possibility to use the linear system of equation (11) for the studying of geo- 
metrical properties of the basic space. 

In particular the invariants of k x k matrix-function 

2 ds 4 

under change of the coordinates ^/f. can be of used for that. 

The first applications the notion of extended spaces for the studying of nonlinear second order 
differential equations connected with nonlinear dynamical systems have been considered in the works 
of author [4, 5, 6]. 

Here we consider some properties of the space defined by the Pfaff equations connected with a 
nonlinear dynamical systems. 



3 The Lorenz dynamical system 

The equations of Lorenz dynamical system are 

dx dy dz , 

— = a{y - x), — =rx-y- zx, — = xy - bz, (12) 

ds ds ds 

where a, b, r are the parameters. 

These equations describe the behaviour of the flow lines along the vector field 

N = [a(y — x), rx — y — zx, xy — bz] 

depending from the parameters a,b,r. 

The properties of non holonomic variety L 2 in this case are determined by the following Pfaff equation 

. .dx dy dz 

vyy - x )-r + \ rx - y - zx )-r + \ X V - bz)— = 0. (13) 
ds ds ds 



The object of holonomicity for the Lorenz vector field is 

(^N, rotN^j — a xy — 2 a x 2 + y 2 — bzr + bz 2 + bza. 



The properties of asymptotic lines of the corresponding variety M 2 are defined by the system of 
equations 

Ts z{s) ) + (lfy P ^ y > z) + lk Q{x ' y ' z) ) (i x{s) ) Ts y{s)+ 

+ (j- g P{x,y,z) + ±R{x,v,z)) ±z{s)+ 



(15) 
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which take the form 



-(^»(-)) +(r + a-z) (±x(8)) ±V{*)-* (^Wl 



d . ,\ d . . , / d 



^i^J^-H^' =0 ' (16) 



(-62 + xy) ^-z(s) + (rx - y - zx) -^-y(s) + (ay - ax) 4~ x ( s ) = °- ( 17 ) 
as as as 



We find from these equations the 



d 1 

i-A s ) = — 



ds a (—y + x) 

4-y( s n rx + (4-y( s n v + (4-A s n zx - (4-A s )) x v + ( 4~A s ) ) bz \^ [ ' S] 



and 
where 



ds I \ds J \ds J \ds J \ds 



A(x(s)) 2 + Bx(s) + C = 0, (19) 



d / N \ 2 ( d , \ ( d , . \ , . , / d 



A=(-az(s) + ar-a)^-y(s)) + ^-y(s) j a ^-z(s) j y(s) - b ^-z(s) j a, 

B = (- (y(s)) 3 + 2bay( S )+y( S )bz(s)) (^*(*)) + 
+ (rbz(s)+2z(s) (y(s)) 2 -abz(s)-(z(s)) 2 b-2r(y(s)) 2 -a (y(s)) 2 + (y(s)) 2 ) x 

x (jis y< ^ s ^) ~<is z ^ + (^ z ^ y ^ +az ^y^^ r2 y^~ ary ^ +2rz ^ y ^ 

-iuy(s)+ry(s)-(z(s)) 2 y(s)) (j^y( s )^J , 

C = (-ba {y{s) f b 2 {z{s) f + (y(s)) 2 bz(s)) (^*{s) 
+ (- (z(s)) 2 by(s) + rbz(s)y(s) - 2 y(s)bz(s) + {y(s)f + a bz(s)y(a)j x 

x (^ y(s) ) Ts z{s) + ( r {v(s))2 {v(s))2 ~ z{s) (v(s))2 ) (^ y(s) ) 2 • 

After differentiating the equation (19) on the variable s and taking into account the expression (18) for 
dx { s ) 

we get the relation 



+ 



ds 



E (x(s)f + F (x(s)f + Hx(s) +K = (20) 
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with some functions E, F, H, K which does not contains the variable x(s). 
The resultant of the equations (19) and (20) with regard of the variable 

d 

d , , d~ s y{s) 
— z(s) = 



ds d 

Tz v{z) 



give us a following conditions 



(r - 1 - z) V(z)^y(z) -bz+ {y{z)f = 0, (21) 



d , V „ d 



L[ Tz y{z)\ +M-y(z) + N = 0, (22) 



where the coefficients of equation are 

L 



= -z% 2 +(2ab 2 + b {y(z)f + 2 rb 2 ) z 3 + 
I — r + z \ / 

+ ((-2 rb - 2 a b) (y(z)) 2 + Aab 2 - r 2 b 2 - a 2 b 2 - 2 a 6 2 r) z 2 + 
+ ((y(z)) 4 + {^b + 2arb + r 2 b-Aab) (y(z)) 2 ) z + (1 - r) {y(z)f , 



-1/2 - M - = -2 z\ 2 +(3ab 2 + 3rb 2 + b (y(z)) 2 ) z 2 + 
y(z) (1 — r + z) V / 

+ ((-2 ab-rb-2b) {y{z)f -2ab 2 r- a 2 b 2 + 4ab 2 - r 2 6 2 ) z+ 
+ (y(z)f + (-3<7b + <7rb + rb + v 2 b) {y{z)f , 



N = b 3 z 4 + (-3 {y{z)f b 2 - 2 rb 3 + 2 a 6 3 ) z 3 + 
+ (b {y{z)f + (5 rb 2 + 4 a b 2 ) {y{z)f + o 2 b 3 + r 2 b 3 - 2 o fe 3 r) z 2 + 
+ ({-rb - 2 a b - 3 b) {y{z)f + (-8 ab 2 r-2 r 2 b 2 + 12 a b 2 - 2 a 2 b 2 ) {y{z)) 2 ^j z+ 

+ {y{z)f + (2rb + a 2 b + 2a rb - b - Aab) {y{z)f + {Aab 2 + Aar 2 b 2 -8crb 2 r) {y{z)f 

The solutions of the first order differential equations (21), (22) together with conditions (16)-(17) 
allow us to get the expressions for asymptotic line of the variety L 2 . 
Let us consider some examples. 

The equation (22) has a set of singular solutions y{z). 
One of them determines from the relation 

z 2 b + (-2 cr b - 2 rb) z + y 2 + 2 bra + br 2 - 4 a b + ba 2 = 

which presents the second order curve in the plane (z,y). 

Remark that the equation (22) presents the algebraic curve 

®(y',y,z) = o 

of genus 3 = 1 with respect to variables y' , y in case 1 and genus g = when r = 1. 

According with general theory some of such type equations can be integrated with help of elliptic 
functions or can be brought to integration of the Rikkati equation. 

In both cases the properties of asymptotic lines should be dependent from the parameters of model. 
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4 The Rossler dynamical system 

Differential equations of the Rossler dynamical model are 

dx dy dz 

— — -y - z, -j- = x + ay, — = b + xz - cz, (23) 
ds ds ds 

where a, 6, c are the parameters. 

For corresponding non holonomic variety V 2 we have a following Pfaff equation 

{-y - 2)4" + {x + ay)^- + (b + xz- cz)% = 0. (24) 
ds ds ds 

The object of holonomicity for the Rossler vector field is 

(jV, rotN^j — —x + xz — ay — ayz + 2 b — 2 cz. 
In this case the properties of the system (16)-(17) for asymptotic lines are determined by the equation 

A^y(z) + B = (25) 



where 



i v{z) 



+ 



A = (—y(z)az + az — az 2 + y(z)a) ^- 
(-2 az 3 - 2 y(z)az - 2 az 2 y{z) - 2 (y(z)f a) ^y(*)- 



-ay(z)z 3 — 2 fez — z 3 c + cz 2 + y(z)c + bz 2 — cy(z)z — a (y(z)) z+ 
+ (y(z) f a + b + az 2 y{z) 



and 



B = 2 (z + y{z)) [\\j- z y{ z )j WW + { ) c " 

Tz y{z) ) a - az {Tz v{z] 

This is equation of the second range at the condition a^O. 
In the case a = it takes a form of the first range equation 

d 2 d 

(2 bz-y(z) c+z 3 c-cz 2 -b+cy (z) z-bz 2 ) — y (z) + 2cz—y (z)+ 

+2 y (z) c^-y (z) + 2 fez + 2 y (z) b = 0. 

At the conditions a = 0, 6 = 0, c ^ the properties of asymptotic lines of corresponding non 
holonomic variety arc dependent from the solutions of the equation 



d^ 1 d d 

(z 3 - z 2 - y (z) + y (z) z) —y (z) + 2 z—y (z) + 2y (z) —y (z) = 0. 
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5 Quadratic systems 

Here we consider the properties of the Pfaff equations connected with a polynomial differential systems 
in R 2 defined by 

dx dy , . ,„„. 

Ts = P (x,y), ^ = «(*,»), (26) 

where p(x, y) and g(x, y) are polinomials of degree 2. 

The system (26) takes the form of the Pfaff equation after extension on a projective plane 

(x Q (x, y,z)-yP (x, y,z))dz - zQ (x, y,z)dx + zP (x, y, z) dy = 0, (27) 

where 

P(x,y,z), Q(x,y,z), R(x,y,z) 

are the homogeneous functions constructed from the functions p(x,y) and q(x,y). 
As example for the system 

^ = kx + ly + ax 2 + bxy + cy 2 , ^ == mx + ny + ex 2 + fxy + hy 2 , 
ds ds 

with a ten parameters one get a Pfaff equation after a projective extension 

P(x, y, z)dx + Q(x, y, z)dy + R(x, y, z)dz = 0, 

or 

((—ny — mx) z 2 + (—hy 2 — fxy — ex 2 ) z) dx + ((ly + kx) z 2 + 
+ (cy 2 + bxy + as 2 ) z) dy + ((—ly 2 + (n — k) xy + mi 2 ) z — cy + 

+ (h-b) xy 2 + (f - a) x 2 y + ex 3 ) dz = 0. (28) 
This equation corresponds the system 

dx dy dz . 

— = P(x,y,z), — = Q{x,y,z), —=R(x,y,z), (29) 

where 

P(x, y, z) = (—ny — mx) z 2 + (—hy 2 — fxy — ex 2 ) z, 

Q(x, y, z) = (ly + kx) z 2 + (cy 2 + bxy + ax 2 ) z, 

R(x, y, z) — (—ly 2 + (n — k) xy + mx 2 ) z — cy 3 + (h — b) xy 2 + (f — a) x 2 y + ex 3 

It is important to note that the vector field N = (P, Q, R) connected with a system (32) is holonomic 
due the condition 

(N,rotN) = 0. 

Corresponding Pfaff (31) equation is integrable and determines the family of developing surfaces 
U(x, y,z) = C in the i? 3 -space. 

The investigation of the asymptotic lines of the surfaces which corresponds the equation (31) may 
be useful in applications. 

Let us consider some of examples. 

The system of equations 

^ = 5x + 6x 2 + 4(1 + n)xy + [iy 2 \ ^ == x + 2y + Axy + (2 + 3^)y 2 , 
ds ds 

with (-71 + 17 v / (17)/32 < u < posseses the invariant algebraic curve 

x 2 + x 3 + x 2 y + 2[ixy 2 + 2[ixy 3 + mu 2 y 4 = 0. 



6 CUBIC SYSTEMS 



10 



The conditions of compatibility of equations for asymptotic lines of the variety defined by the 
corresponding Pfaff equation lead to a following conditions on the functions 

(3 [? + 4 fi + 2) (y (s)) 2 + (9 n + 6) x (s) y (s) + 6 (x (s)) 2 = 0, 

(y (s)) 2 n+(-2n-l)x (s) y (s) - 3 (x (s)) 2 = 0. 
The values of parameter \x determined by the conditions 

(3/i+10)(M + 4)(/i-2) = 0, /i = 

are special. 

Next example is the system with at least a four limit cycles. 

dx o 2 dy o r 

— = kx - y - lOx + bxy + y , — == x + x + fxy, 
ds ds 

The studying of asymptotic lines for this system give a following conditions on parameters 
10 k 2 + 11 fk 2 + fk 2 - 7 k + 2 fk - 9 b + k 3 - 80 - 81 / + b 2 k + 2 bk 2 + 9 bk + 9 fbk = 
and functions 

x 3 + 10 y (x) x 2 - (y (x)) 3 + fx 2 y (x) - bx (y (x)f = 0, 
/ (y (x)) 2 — fy (x) kx + (y (x)) 2 +y(x)bx + y (x) x — kx 2 — 10 x 2 = 0. 
Remark that these equations equivalent the equations of direct lines. 

6 Cubic systems 

By analogy can be investigated the properties of the asymptotic lines of the surfaces connected with 
the system 

dx dy 

— =p(x,y), —=q(x,y), (30) 

where p(x, y) and q(x, y) are polinomials of degree 3. 
Let us consider the examples. 
The system 

dx dy 2 2 

Tg =y, - = -x-xy + muy (31) 

connects with the Van der Pol equation. 

After extension on the projective plane we get an integrable Pfaff equation 

(— x 2 z 2 — x 3 y + x[i 2 yz 2 — y 2 z 2 ) dz + (z 3 x + zx 2 y — z 3 [i 2 y) dx + yz 3 dy = 0. 

The equations for the asymptotic lines of corresponding surface give the conditions 

x 2 + (y (x)f -x/j, 2 y(x) = 

and 

— x 2 + 2 (y (x)) 2 + 4x^i 2 y (x) = 0. 
From the first condition we find 

y(x) = (l/2 1? + (-)l/2 Vm 4 - 4) x, 

and the second gives us 



V(x)= (-M 2 + (-)l/2 V / 4^ 4 + 2 



x. 
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For the system 

^ = -y + ax(x 2 + y 2 - 1), ^ = x + by(x 2 + y 2 - I) (32) 

the point (0,0 is node at the condition ab > — 1 and (a — b) 2 >= 4 and it has the limit cycle around 
this point. 

After extension we get an integrable Pfaff equation 

[x 2 z 2 + byx 3 + xby 3 — xbyz 2 + y 2 z 2 — yax 3 — axy 3 + yaxz 2 ) dz+ 

+ (—z 3 y + zax 3 + zaxy 2 — z 3 ax) dy + (—z 3 x — zbyx 2 — zby 3 + z 3 by) dx = 

with a developing surfaces as general integral. 

The conditions on parameters for existence of asymptotic lines are 

a6+l = 0, -a + b = 0. 

For the corresponding functions one get the equations of direct lines 

xby (x) — x 2 — y (x) ax — (y (x)) 2 = 0, 

or 

y(x) = xby (x) - x 2 - y (x) ax - (y (x) f , 

and 

(30 ab + 8 a 3 b - 25 a 2 - 9 b 2 - 12 a 4 ) x 6 + 
+ (36 b 3 + 76 ba 2 - 108 ab 2 - 4 a 3 - 32 a 3 b 2 + 32 ba 4 ) yx 5 + 
+ (69 b 2 + 53 a 2 + 36 a 4 - 64 a 3 b - 134 ab + 24 ab 3 - 8 a 2 b 2 ) y 2 x 4 + 

+ (-24 b 3 - 8 ba 2 + 64 a 3 b 2 + 8ab 2 + 24 a 3 - 64 a 2 b 3 ) y 3 x 3 + 
+ (-8 a 2 b 2 - 134 ab + 24 a 3 b + 36 b 4 - 64 ab 3 + 69 a 2 + 53 b 2 ) y 4 x 2 + 
+ (-32 b 4 a - 36 a 3 - 76 a6 2 + 32 a 2 b 3 + Ab 3 + 108 ba 2 ) y 5 x+ 
+ (-25 fe 2 - 9 a 2 + 30 ab + 8 ab 3 - 12 b 4 ) y e = 
with a very complicate relations between the parameters a, b. 

7 Quatric systems 

The system 



^ = -Ay + yx 2 - x 4 , $=as-i 3 (33) 
as as 



has a center and the limit cycle at the conditions a = 2 A + yl 2 ,0<A<5.10 5 . 
After extension we get an integrable Pfaff equation 

(z 2 x 3 — z 4 ax) ^—x (s) + (— x 4 z — z 4 Ay + z 2 x 2 y) ^—y (s) + 
ds ds 



+ (ax 2 z 3 - x 4 z + Ay 2 z 3 - zx 2 y 2 + yx 4 ) —z (s) = 

ds 

with a developing surfaces as general integral. 

For the functions one get the equations of direct lines 

x 10 + (3 a - 9 A) z 2 x 8 + (24 A 2 - 3 Aa) z 4 x e + (A 2 - 36 A 3 - 2 Aa + a 2 ) z 6 x 4 + 

(6 A 3 + 36 A 3 a - 12 A 2 a + 6 a 2 A) z 8 x 2 + (9 a 2 A 2 - 18 A 3 a + 9 A 4 ) z 10 = 
and the family of conies 

(z( S )fa-(x( S )) 2 = Q 
(x(s)) 2 -A(z( S )) 2 =0. 



8 GEODESICS OF THE FIRST KIND 12 



8 Geodesies of the first kind 

A geodesies of the first kind on non holonomic variety are denned by the system of equations (5)-(7) . 

Let us consider the solutions of this system of equations for the variety V 2 defined by the vector 
field 

F = [ayz, bxz, cxy], 

where a, b, c are parameters. 

In this case the condition N = holds and variety is holonomic. 
The system of equations for geodesies is 

d2 ayzx (c + b) (j^yi 3 )) T" z ( s ) a V 2z ( a + c ) ( T" 3 ^)) ~d~s~ Z ^ 
x(s) H - - I - - h 

ds 2 a 2 y 2 z 2 + b 2 z 2 x 2 + c 2 x 2 y 2 a 2 y 2 z 2 + b 2 z 2 x 2 + c 2 x 2 y 2 



ayz {a + b)^-x(s)j ^g(«) _ o 
a 2 y 2 z 2 + b 2 z 2 x 2 + c 2 x 2 y 2 

d 2 bx 2 z (c + b) ^-z(s) bxzy (a + c) ( ±x(s) j ^z(s) 

y( s \ -| (Ml | 2 L |_ 

ds 2 a 2 y 2 z 2 + b 2 z 2 x 2 + c 2 x 2 y 2 a 2 y 2 z 2 + b 2 z 2 x 2 + c 2 x 2 y 2 



bxz 2 (a + b)(f s x(s)) —y{s) 
a 2 y 2 z 2 + b 2 z 2 x 2 + c 2 x 2 y 2 



d 

Vis) 

0, 



d 2 „ cxy(c + b) ^-y(s) j -z(s) cxy (a + c) ^-x(s) j -z(s) 
ds 2 a 2 y 2 z 2 + b 2 z 2 x 2 + c 2 x 2 y 2 a 2 y 2 z 2 + b 2 z 2 x 2 + c 2 x 2 y 2 



cxyz(a + b) [^ X ( S )J ^v( s ) 
a 2 y 2 z 2 + b 2 z 2 x 2 + c 2 x 2 y 2 



This system of equations has an integral 



ay(s)z(s)-^-x(s) + bx(s)z(s)-^-y(s) + cx(s)y(s)-^-z(s) = 
ds ds ds 

and can be integrated. 

In fact after substitution of the expression j^z(s) into the above equations one get the algebraic 
relations with respect the variable z(s) which are compatible at the conditions 

d ( d \ ^ 

(ab + b 2 ) {z{y)f + 2 bz(y)cy—z(y) + (ay 2 c + c 2 y 2 ) i^ z (y)J = 

and lead to the solution 

z(y) = y K _Cl 

where 



cb yt—cab (a + c + b) 
c 2 + ac c 2 + ac 

Finally we present the expression for the Chcrn-Simons invariant of affine connection defined by the 
equations of geodesies (5-7) 

CS = J e^ k {T%Tl p . 3 + 2 -T%Y) r Tl p )dxdydz. (34) 
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In the case of the Lorenz system of equations it has been obtained with the help of a six dimensional 
Riemann extension of corresponding space and has the form 



L = (26+2 - 2(j)x 2 y 2 + (3 a 2 + 4ar - 4r6 - 2ba - 4r) zx 2 + (26 + 2 - 2 a) z 2 x 2 + 

+ (-3 ra 2 + 4 er 2 6 - 5 a 3 + 2 6cr r - 2 a r 2 +4 a 2 + 2 br 2 + 2 r 2 ) x 2 + 
+ (-2 rb - 4 r + 9 a 3 + 2 ra 2 ) yx + (-2 a + A - 2 a 2 - 2 a r + 2ba - Ab 2 ) zyx- 
- ((bar + 2 a 2 + a 2 b - a z 2 - 2a r - a r 2 ) y - ay 3 ) x + (-2ab 2 + 2b 3 - 2b 2 ) z 2 + 
+ (-4 a 3 - a 2 b + 2 - 2 b - a 2 + a r + (-a + ba) z - ba r - 2 a) y 2 + 
+ (-2 a b 2 + 2 rb 2 + 2 b 2 ra - 2 b 2 a 2 ) z, 
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where 



and 



M = (x 2 + b 2 ) z 2 + ((-2 6 + 2) xy - 2 rx 2 ) z+ 
+ (a 2 + 1 + x 2 ) y 2 + (-2 a 2 - 2 r) xy + (a 2 + r 2 ) x 2 . 



